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Abstract 

We use a light-cone sum rule (LCSR) analysis of the branching ratios of radiative meson decays to 
contrain the value of the magnetic susceptibility of the quark condensate x(/-0- F° r me fi rst time, 
we use a complete set of three-particle distribution amplitudes that enables us to give a consistent 
prediction for the branching ratios. Furthermore we will make use of a very recent update of several 
non-perturbative parameters. Our final result for = 1 GeV) = 2.85 ± 0.5 GeV~ 2 (assuming 
asymptotic wave functions) agrees with the currently used value of 3.15 ± 0.3 GeV -2 . 



1 Introduction 



The term "magnetic susceptibility of the quark condensate", was coined as early as 1984 in a 
fundamental work by Ioffe and Smilga HI on the magnetic properties of the QCD vacuum. The de- 
notation is due to a very intuitive picture: the value of the non-perturbative parameter x determines 
the strength of the reaction of the vacuum with respect to an applied external electromagnetic field, 
analogous to what is understood by magnetic susceptibility in solid states physics. 

The magnetic susceptibilityLl is an important input parameter in various calculations of pro- 
cesses involving real photons, such as dijet photo-production [2], radiative decays [3] or the anoma- 
lous magnetic moment of muons [4]. It is known that x is i n f act surprisingly big compared to what 
one might assume as its "natural" scale, hence, it often provides the dominant contribution to the 
calculation, see (51 for examples. Therefore, while the determination of the magnetic susceptibility 
may provide additional insights into properties of the QCD vacuum, a task worthy in its own right, 
a precise knowledge of the value of x(aO is important for its application of QCD calculations. 

Thus, various theoretical approaches to determine the magnetic susceptibility have been sug- 
gested. To the best of our knowledge, the first work regarding this topic was published in 1983 
[6|| and made use of the so-called vector-dominance approximation. Others include consistency 
conditions of QCD sum rule predictions [1] and local duality approaches [H. The value of 
XifJt = lGeV) = 3.15 ± 0.3 GeV -2 , which is currently used in most QCD sum rule [9] calcula- 
tions was obtained by a combination of a local duality approach and a QCD sum rule for \ OH- 
However, in [11] the possibility of a larger value is considered using a rather elegant argument 
similar to the famous Gell-Mann-Oakes-Renner (GMOR) relation. An approach based on the 
instanton liquid model can be found in lfT2l[T3l[T4l . 

In this work a different way to determine the value of x is employed. We fit a %-dependent 
sum rule to an experimental value. We will argue (in Sec© that radiative heavy meson decays 
and especially the D* — > D'j decay and the corresponding branching ratio are suited for such a 
procedure. Heavy meson decays have been studied numerously and successfully with sum-rule- 
based approaches, see e.g. [15]. In lfT6ll the different D* — > D^f and B* — > B^f decays have 
been examined using a light-cone sum rule approach based on the background field method, that 
was advocated in IfTTl . While we follow the general spirit of this calculation, we use a more 
recent complete set of photon distribution amplitudes up to twist 4 ifTOl . including some small 
contributions that where neglected in lfT6ll . We will also make use of a recently published update 
on some non-perturbative parameters iTTHl . As a detailed calculation of the various D* — > Dn 
decays was presented in Q~9] and the D* decays only into Dn and D'j, we are able to determine 
the ^-dependence of the corresponding branching ratios up to and including twist 4. 

The paper is organized as follows. In Chapter [2] we will calculate the D* + — > D* +r y decay 
constant g D * Dj using LCSRs. The next Chapter deals with the branching ratios and how to extract 
the value of X- Chapter |4] contains the conclusions. 

*As we only deal with one kind of "magnetic susceptibility", we will usually drop the reference to the quark 
condensate. 



1 



2 LCSRs for the D* L> 7 decay 

As already mentioned in the introduction, if one wants to find \ by fitting a sum rule to experiment, 
it is important to select a process the exhibits a strong dependence on x- This unfortunately rules 
out some well known classical problems like the magnetic moments of nucleons, as sum rule 
calculations [5] show large contributions due to higher twist terms, most notably twist 4, and only 
a moderate dependence on X- A suitable process should thus fulfill several conditions: 

• in order to archive a better suppression of higher twist terms it is advantageous to consider 
processes that involve a propagating heavy quark, which will induce a — — suppression of 
twist-4 contributions 

• existing experimental data should be accurate enough to support a fit 

• the relevant sum rule should not contain any theoretical difficulties, e.g. sum rules for tran- 
sition between states of different mass are known to require an adjustment of the Borel 
parameters to allow for the mass gap, which would possibly entail additional uncertainties. 

This does not single out any particular process, however our choice, the D* — > decay, seems to 
fulfill the criteria rather well \ It is also advantageous that the decay D* — > Dtt has been analyzed 
using LCSRs lfT9l , thus allowing us to consider sum rules for the branching ratios, which have 
been measured rather precisely If20l . 

In the remainder of this Chapter we will examine the D* — > D^f transition. In Ref. |fT6l this 
was already calculated using LCSRs in conjunction with photon distribution amplitudes. Our 
calculation will use, for the first time, a complete set of twist 4 DAs and updated parameters. 
As the only difference between the D*° and the D* + is the charge of the light quark, we will focus 
on the case of the D + §. 



2.1 Definitions 

The transition matrix element 

(D*(p, A) | j^(0) \D{p + q)) = e^p u 6 {X)a q p g D * Dl (Q 2 ) (1) 

can be parametrized by introducing the transition amplitude go* o 7 . Here e> x >P is the 4-polarization 
vector of the D*. The decay width only depends on 5 f D*D 7 (0), therefore, we only have to take 
real photons into account. Henceforth, ej* is the four-polarization vector of the emitted photon, 

q ■ eW = 0. 

The decay as shown in Fig Q] can conveniently be described by the correlation function 

IL^ip,q) = i 2 jd'x J^ye f ^(0\T{r^(x)j v ( V )li D (0)}\0) . (2) 

tin fact, at first glance the decay B* — * £>7 is even more appealing. Unfortunately the B* is too heavy to be 
produced in B-factories running at the T(4s) resonance. Furthermore, due to kinematics the process B* — > Bit is 
only possible as a virtual subprocess, therefore the branching ratio is not a good observable for our purpose. 

*To avoid unnecessary indices, we will usually abstain from using the + . 
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Figure 1: Radiative decays of a D* with momentum p into a D with momentum p + q. 

Here 

= eadj^d + e c C7 M c (3) 
is the electromagnetic current, with = — 1/3 and e c = 2/3 being the quark charges. The currents 

ma 

v ^(x) =v a (x)Yc a (x) (4) 

77^(x)=#(x)i 75 c 6 (a;) (5) 

are used to generate a state with the quantum numbers of a D* or a D meson, respectively, a, b are 
color indices. The coupling constants fr>* and f D of the currents © and © to the corresponding 
states are defined as [fTTl 

<0| 77^(0) \D*(\,p))=m D *f D *e^ (6) 

9 

(0\ VD (0)\D(p))=^f D . (7) 

m c 

Hereafter, m D * and m D is the mass of the D* and the D. m c is the mass of the charm quark. 

It is advantageous to rewrite the correlator by introducing an electromagnetic background field 
of a plane wave 

F,u = i {el x k» - efq v ) e**> . (8) 
This allows us to write the following correlation function 

K U (P, q) el X) = i Jd i xe^(0\T{r ] ^(x)rJ D (0)}\0) F . (9) 
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Here the subscript F indicates that the vacuum expectation value has to be evaluated in the back- 
ground field F^ LV . The correlation function in Eq.© can then be reproduced by expanding Eq.® 
in powers of the background field and taking only the terms linear in F^ v corresponds to the single 
photon emission. A detailed analysis of the general procedure can be found in Ref. lfTOl . while an 
excellent review on the background field method can be found in [|2TT|. 



2.2 Sum rule for g D * Dl 

Following the general strategy of QCD sum rules, the correlation function Eq.© has to be evalu- 
ated in two different kinematic regions. On the one hand, Eq.© will be dominated by the decay 
D* — > D'j, if p 2 ps m 2 D * and (p + q) 2 ps m 2 D . Then again, in the kinematic region where p 2 <^ 
and (p + q) 2 <C an expansion in terms of photon distribution amplitudes of increasing twist is 
valid. 

On the hadronic level Eq.© can then by written as 

n"" (v a) e (A) - fp*fD m D m D*9D*D~t m&p (x )p (m 
11 P [p, q) e v - ^ ^ _ ^ ^ _ p a qpe +... (1U) 

Where we made use of Eqs. dl 16171) and introduced the abbreviations Px—p and p 2 = p + q. The 
dots represent contributions from excited states and the continuum. 

To calculate the correlation function in the deep Euclidean regime is is necessary to insert the 
explicit expressions for the currents rf D * and r/ D into Eq.©, yielding 

rP"(p,g)4 A) = % J d 4 xe ipx (0\T \d a {x)Yc a (x)c b (0)i-f 5 d\0)} \0) F . (11) 

As both massive and massless quarks propagate in the presents of gluonic and electromagnetic 
backgroundfields the following expressions^ have to be used after employing Wick's theorem 11221 : 



d(x)d(0) = - / du {ufo a p + ua a p$} G ap [ 

Z7T X l07T X In 



ux) 



1Qtt 2 x 2 J 

and 

c(x)c(0) = / — — e 



I du {u$a a p + uo a p$} F a/3 (ux) + . . . , (12) 
Jo 



(2tt)H m 2 c -k 



2 



d A k 

(2?r)V 



-ikx 



du 



+ m c 



2 m 2 



k 2 ) 



; G^{ux)a^ + 



ux. 



m 2 . — k 2 



G» v {ux) lv 



§We tacitly assume Fock-Schwinger Gauge, which allows us to omit the path-ordered exponents. 
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diagram a 



diagram b 





diagram c 



diagram d 



Figure 2: Diagrams up to twist 4. The wiggled and the curled lines represent the coupling to the 
electromagnetic and gluonic background fields. The double lines denote the heavy quark and the 
single lines the light quark. Photon distribution amlitudes are indicated by crosses. 



- ie c 



d 4 k 



-ikx 



du 







+ m c 



-2\2 



ux, 



m 



+ ... 



(13) 



where we used the common abbreviation u = 1 — u. F^ LU (x) is the electromagnetic field strength 
tensor defined in Eq.® and G^ v = G A v t A is the gluon field strength tensor. The dots represent 
terms that will not give rise to terms of twist 4 or lower. 

Up to twist 4 there are only four relevant Feynman diagrams as depicted in Figj2] and their 
counterparts with exchanged heavy and light quarks. However, as the charm quark condensate 
vanishes identically, only the perturbative "mirror diagram" has to be considered. The calculation 
is then straight forward, although quite cumbersome, so that we will just give the final results for 
the different diagrams. 

The sum of the two perturbative diagrams is given by 



rpfl 



An 2 



dt 



te r m r 



+ t-. 



te d m c 



tm 2 , — tt{up\ + up 2 ,) tm 2 — tt{up\ + up 2 . 



/j.va/3 (A) 



For the remaining diagrams one finds: 



Tbip, q) 



£d{qq) I du 
o 



m 2 — {up\ + upl) 



(14) 



A(u) 



+ 2 



m 2 — (upl + upl)) 2 (m 2 — (upi + upl)) 3 



+ e d m c 



du 



2(m 2 — (upl + upf))' 



^ A (u) 



Pvqa&p 



(A) 



(15) 



5 



t c(p, q) = e c (qq) 



Uo ^ {ml - (p + a u q) 2 ) 2 



(5 7 (o) - (1 + 2 M )T 7 4 (a)) + 2P [T 7 4 («)j 



(16) 



1 

V -{m 2 c -(p + a u q) 2 ) 2 



S{a) + (1 + 2u)(T 3 (o) - T 4 (a) - 5(a)) 



- 3Ti(a) + 3T 2 (a) + 2P [Ti(a) - T 2 (a) - T 3 (a) + T 4 (a)] 



fivafl 



P»q a e 



(A) 

(17) 



Here J Da = da q da-g J da g 5(1 — a q — a-q — a g ). The functions <p and ip A , which 
are of twist 2 and twist 3, respectively, and Tj, 5, 5, <S 7 , Tj and A, which have twist 4, are 
defined in Appendix [A] In order to simplify the expressions, several rather lengthy terms have 
been abbreviated via the function F[x] defined in Appendix |B] 

The results for the diagrams Figj2^,b reproduce those given in lfT6l . if one takes into account 
the different nomenclature for the DAs and the different sign convention for the electromagnetic 
current. The contribution being rather small was not taken into account in IfToll . 

In order to obtain a sum rule for go*D-yi it is necessary to equate the hadronic representation 
Eq. (flOl) and the twist-expansion, Eqs. (ll4ll5ll6ll7l) . of the correlation function. This requires to 
control the contributions of excited states and of the continuum, that can be approximately evalu- 
ated using quark-hadron duality. The standard procedure involves a Borel transformation, which 
suppresses these contributions and a subsequent continuum subtraction. 

As the momentum of the photon q does not vanish in the real photon case q 2 = 0, we can treat p 1 
and p 2 as independent, which allows a so-called double Borel transformation. This transformation 
can be carried out with the help of the formulas 



B M? B 



Ml 



T(a) 



up 2 2 ) 



t 2 ~ a 6 u - 



M 2 



M 2 + Mf 



(18) 



and 



B Ml B 



Ml 



(m\ - pj) [ml - p 2 2 ) 



-m\/Ml-m%/Ml 



(19) 



Here Mf is the Borel parameter corresponding to p 2 and t = MfM|/(i\if + Mf). As our final 
aim is a prediction for the branching ratio with help of the results from [19], the method of how 
the continuum is subtracted must match the one employed therein. We will not go into the details 
of the procedure but instead would like to refer to the corresponding sections in [|T9ll for a more 
detailed discussion. Furthermore, it should be noted that the non-trivial continuum subtraction of 
the perturbative contribution, Eq. (fT4l) . is elaborated in |fT6l[T9l . 



6 



This leads to the following sum rule for g D * Dy : 



e * 



9d*+d+ 



7 



f D *fDm 2 D m D * 



m c N c 
Att 2 



So 



ds 



m;\ I s 

+e c ln — 
s J \m 2 



-s /t 



.so. 



(-d 



^^(f) A (v) + e d tx (qq) y?0) 



with 



A(u) ,_ v 
e d^~ (qq) 



1 + 



m; 



M 



, = 2m 2 D ,m 2 D 
m 2 D + m 2 D * 



G 



(20) 



(21) 



2s 2 n *S 



2 

D*°D 



(22) 



M| 



(23) 



M| + M ? ' 

S£> and S£>* are the continuum thresholds for the D* and D, respectively. The functions X F and X G 
correspond to lengthy contributions of twist 4. The full expressions can be found in Appendix E 
The decay width T(D* + -> D+7) is then given by the formula AH 



T(D* + -> L> + 7 ) = 

The corresponding expressions for the r(D*° 
<-> e u in Eqs. (|2QI24l) . 



m 



(24) 



384tt 2 

-> -D°7) decays can easily be found by replacing 



3 The branching ratios and the determination of x(/i = lGeV) 



The branching ratio of the D* + — > D +r y decay is defined as 



B (£>*+ -> D+ 7 ) 



r(z> 



D+7) + r(D*+ -> D+7i°) + r(X>*+ -> d°tt^ 



(25) 



In order to determine B (D* + — > £> + 7) one needs the expressions for the decay widths of the 
/}* — > D7r decays. Up to twist 4 accuracy these have been determined in |[T9l 0, which allows us 



more recent evaluation Il24l also includes a s -corrections, for consistency we do not take these into account. 
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Figure 3: The left panel shows the plot of the branching ratio B(D* + — > D +r y) versus the Borel 
parameter t for different values for the product % (//)<£> (1/2). The solid line represents the best fit 
to experiment and uses x(A t )^(l/2) ~ 5.1 GeV -2 . The dashed line corresponds to x(//)</?(l/2) ~ 
6.5 GeV" 2 , the dashed-dotted to x^MV 2 ) ~ 4.5 GeV" 2 . The right panel shows the plot of 
the branching ratio B(D*° -> D ^). The best fit (solid line) yields x(^)y(l/2) « 3.5 GeV" 2 . 
The dashed and dashed-dotted lines correspond to x( / u)(/?(l/2) ps 4.8 GeV" 2 and x(//)</?(l/2) 
6.5 GeV" 2 , respectively. 



to construct a sum rule for the branching ratio using Eq.([20T) and equations (44), (52) and (53) from 
Ref.[19|. A single sum rule for B (D* + — > D +r y) has the distinct advantage that the dependence on 
several parameters that effect the individual sum rules for the widths will be significantly reduced. 
Most notably the dependence on the coupling constants fu and fu* , which are known only up to 
at best 10%, will vanish completely. It is, in principle, not possible to distentangle % from ip(v). 
Hence, we can only make a prediction for the product of the two. 

As the mass difference m D * — m D is very small (of the order 0.07m£)*) and can therefore be 
neglected, it is possible to use symmetric Borel parameters M 2 = Mf , which corresponds to v = 
1/2. In the following analysis we will use the continuum threshold So = 6 GeV 2 as determined 
from the sum rules for the couplings f D and f D * [fT9l , a charm quark mass m c = 1.3 ± 0.1 GeV 
and the Borel window 2 GeV 2 < t < 4 GeV 2 . The experimental value for the branching ratio [|25l 

B (D* + -> D+ 7 ) = 1.6 ± 0.4% (26) 

will be an additional input parameter. In Fig. [3] the plot of our sum rule for the branching ratio 
is shown using various choices for x(v — 1.3GeV)<^(l/2). It can be seen, that branching ratio is 
indeed very sensitive to the value of this product. The best fit is achieved with 

+(0.4+0.3) GeV - 2 _ 



[x^(l/2)](^ = 1.3GeV 2 ) = 5.i; 



(27) 



-(0.7+0.3) 

The first given error is due to theoretical uncertainties, the second error stems from the experimen 
tal bounds. 



The branching ratio for the D ' — > D u 7 decay 



B (D*° -> D° 7 ) 



Y(D*° -> D° 7 ) + T(D*° -> D°tt) 
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can be determined analogously and the corresponding plot can be found in Fig|3] The sum rule has 
a weaker dependence on the fit parameter compared to the previous one. It should further be noted 
that the experimental uncertainties are significantly smaller. A value of 

fop(V2)](/< = 1-SGeV 2 ) = 3.5+;™;*] GeV' 2 (28) 

produces the best agreement with experiment 11251 

B (D*° -> D°j) = 38.1 ± 2.9% . (29) 

The identification of the errors follows Eq.(l27l). It should be noted that the result in Eq. (|28l shows a 
rather strong dependence on the charm quark mass, which is one reason for the rather large relative 
errors. However, this is not the case for Eq.d27l) as the sum rule for B (D* + — ► D + 'j) is almost 
insensitive to the c-quark mass. 

In order get an expression for x alone one has to introduce an explicit value for ip(l/2). The 
natural choice is to use the so called asymptotic form (p(u) = Quu, see discussion in [[Toll . Taking 
into account the scale dependence of the product x<P, one gets 

X (fi = lGeV 2 ) = 3.5±g;| GeV" 2 (30) 

from the results for D* + — > D +i y and 

X(p = lGeV 2 ) = 2At° 9 7 GeV~ 2 , (31) 

form D*° -> D° 7 . 

For completeness, we also determine the decay widths for the different radiative D-decays. In 
this case on has to use a fixed value for = 1.3GeV)<^(l/2), we use 4.73 ± 0A5GeV~ 2 , which 
corresponds to <p(l/2) = 3/2 and xO = lGeV 2 ) = 3.15 ± 0.3GeV" 2 . In this case, the coupling 
constants fjj and fo* are relevant and we will use the sum rules given in Eqs. (46,47) in [|T9l . The 
results are shown in FigJH 

T(D*° -> D°j) = 20 ± 6 keV (32) 



T(D* + -> D+7) = 0.55 ± 0.3 keV (33) 

The large uncertainties are mainly due to the uncertainties in the coupling constants f D and f D * , 
that are not present in the sum rules for the branching ratios. A comparison with the CLEO data 
11201 r(D*+) totaX = 96 ± 26 keV would lead to a branching ratio of w 0.6%, which would be below 
the value given in (|26l) . While the experimental number has yet to be determined more precisely, it 
seems to indicate that QCD sum rules generally underestimate the decay widths in D* — > and 
D* — > Dn, see also [|24l[T6"ll23~1l . However, the LCSR predictions for the branching ratios, that are 
measured more precisely then the widths, typically agree rather well with experiment. 
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t in GeV 2 t in GeV 2 

Figure 4: The sum rules for the radiative D* decay widths using %^(l/2) = 4.65 GeV~ 2 . 

4 Conclusions 

We have calculated the radiative decay constants of the D* using the approach of light-cone sum 
rules. This enabled us to use results from |fT9l to construct a sum rule for the branching ratios 
B(D* + — > .D + 7) and B(D*° — > D 0r j). As this sum rule is dominated by terms of twist 2, that are 
proportional to the magnetic susceptibility of the quark consdensate x> it is possible to determine 
its numerical value by fitting the sum rule to the experimentally determined value of the branching 
ratio. However, the magnetic susceptibility cannot be accessed directly as it is "masked" by an 
additional factor ip(l/2), and only the product of the two can be fitted. 

Taking the naive average of Eq. (|27l) and Eq. (|28~l) . the best agreement with experimental data 
can be achieved by the choice 

[Xf(l/2)](n = 1.3GeV) « 4.3 ± 0.7 GeV" 2 . (34) 

In order to be able to compare our result to other calculations of x, it is necessary to assume an 
explicit value for cp(l/2). The choice of the asymptotic wave function, which has produced good 
agreement with experiment in the past, leads to 

X(ji = IGeV) = 2.85 ± 0.5GeV" 2 . (35) 

This value is below the most recent QCD sum rule result IfTOl and is still affected by rather large 
uncertainties, which are due to the limited precision of both sum rule parameters and the experi- 
mental data. However, this result is very close to the one obtained in [14] using an instanton liquid 
model. Alternatively, assuming the standard value = 1 GeV) = 3.15 GeV -2 the estimate 

(p{l/2, ii=l GeV) « 1.35 ± 0.3 (36) 

can be obtained. At first glance, this seems to indicate that the shape of the wave function may be 
more fiat compared the the asymptotic form, however, our calculation is not yet precise enough to 
support such a statement. 
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A Photon distribution amplitudes 



For completeness we collect the relevant photon distribution amplitudes for the D* — > D^y tran- 
sition according to tfTOl . Note that in ifTOl the photon momentum has the opposite sign and the 
parameterization of the separation of antiquark and quark is different. 
The path-ordered exponents 

[x, y] = Pexp jz Jdt(x- y)^ [e q A^(tx - ty) + gB^(tx - ty)] j 

assure gauge invariance of the matrix elements. It is important, that the electromagnetic potential 
is included in addition to the gluon potentials^, as additional terms to those given in IfTOl will 
occur otherwise. 



A.l Twist-2 and Twist-4 DAs 

The leading-twist DA reads 

(0|g(0) [0,x]a a/3 q(x) \0) F = e q (qq) / du x<f(u)F a p(ux) 

Jo 



16 In 



with 



+ 6 -^p- [ du B(u)af {xpF ap {ux) - x a F Pp {ux)) (37) 
° Jo 



<p( u ) = ^\u) =Qu(l-u) 



(38) 



A(u) = 40w(l - u) (3k - k + + 1) + 8 (C 2 + - 3C 2 ) 



x 



[m(1 - u) (2 + 13u(l - u)) + 2m 3 (10 - 15m + 6m 2 ) hx(u) 



+2(1 - uf (10 - 15(1 - u) + 61 - u 2 ) ln(l 



u 



(39) 



pU 

B(u) = 40 / da (u - a) (l + 3k + ) 
Jo 



(40) 



(0\q(0)e q [0,x] F^(ux)q(x) \0) F = e q (qq) / X>a5 7 (a)F MI/ (a u x 



(41) 



12 



(0\q(0)e q [0,x] a aP F txv (ux)q(x) |0) F = 



e g {qq) 



qx 



[QaQ^e^xp - qp%e\ v - q^e^x p + qpg^e^Xa] 7?(u, qx) 



(0\q(0) [0,ux] gG^iux) [ux,x] q(x) \Q) F = e q (qq) VaS(a)F^ u (a u x 



(0\q(0) [0,ux] i^gG^ (ux) [ux,x] q(x) \0) F = e q (qq) JvaSio^F^otyx) 



(0|g(0) [0,ux] a af3 gG^(ux) [ux,x] q(x) \0) F 



-e q (qq) q a e { ^ v - qpe^g^, - q a ef u g^ + Ti(u, qx) 



e g (qq) 



(A) _L (A) _L (A) _|_ (A) J_ ' 



T 2 (u,qx) 



e q (qq) 



qx 



[(laq^pXv - qpq^e\ a - q a q v e]_ p x „ + g/jg„e* a a;J T 3 (u, qx) 



e q (qq) 



qx 



[q a q»e x ±u xp - qpqpe\_ v - q a q v e\^x p + qpq v e\^x a ] T 4 (u, qx) 



used 



/Da = / da. q \ da q / da g 5(1 — a q — a q — a g ) 
Jo Jo Jo 



9uu 9iiv 



qx 



ei W =9^ 



5(a) = 30a 2 g [(« + k + ) (1 - a g ) + (Ci + Ci + ) (1 - <* g ) (1 - 2a g ) 
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+C 2 (3 (a ? - a q f -a g (l- a g ))] 



(50) 



with 



5(a) = -30a 2 g [(« - k+) (1 - a g ) + (d - Ci + ) (1 - %) (1 - 2a g ) 

+C 2 (3 (a q - a q f - a g (1 - a g ))] (51) 

5 7 (a) = 60«g (a ¥ + a g ) (4-7 (a g + a^)) (52) 

2i(a) = -120 (3( 2 + C 2 + ) (a? - "<?) a q a q a g (54) 

T 2 (a) = 30a 2 g (a q - a,) [(« - «+) + (Ci - Ci + ) (1 - 2a 9 ) + C2 (3 - 4a,)] (55) 

73(a) = -120 (3( 2 - C2" ) ( a q - a q) Oiqa q a g (56) 

T 4 (a) = 30^ ( % - a q ) [(« + k+) + (Ci + Ci + ) (1 - 2a 9 ) + C2 (3 - 4a s )] (57) 



T 4 7 (a) = 60a 2 g (a q - a,) (4-7 (a, + a?)) . (58) 

The abbreviation a represents (a q , a q , a g ). The values of the various constants can be found in 
tabled! 

It should be noted that the matrix element 

(0| q(0)e q [Q,x] a aP F fJLV (ux)q(x) |0) F 

vanishes exactly if one sums up the whole conformal expansion. The expansion itself has, however, 
non-zero coefficients and thus in next-to-leading order in conformal spin the matrix element is 
different from zero. For the same reason the matrix element 

(0|g(0)e, [0,x]F^{ux)q(x)\0) F 

has herein mentioned form and not e q (qq) F^ v (ux) . 



A.2 Twist-3 DAs 



/1 
du^ {V \u)x"F pa (59) 
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A.3 Numerical values for the parameters at the renormalization scale \x = 

lGeV 



A. 


3.15 ± 0.3 GeV" 2 


K 


0.15 


K+ 


-0.05 


Cl 


0.4 


ct 





C2 


0.3 


C 2 + 







-(4 ±2) ■ 10- 3 GeV 2 


«f 


-2.1 ± 1.0 




3.8 ± 1.8 


(9?) 


-(240 ± lOMeV) 3 



Table 1 : Numerical values and uncertainties of the relevant parameters IfTOl [i~8l [T71l 



B The functions P and X F , X { 



G 



In this Section we have gathered the explicit expressions for the three functions P, Xp and Xq that 
appear in Eqs. (|16|17t and (1201 



P[X(a)} =8e c (qq) / dw / da q / dc^ 
Jo Jo Jo 



da-aU 



pqX(a q , a q , 1 - a q - a q ) 



TOO 



m c — (P + ( M — Ua q + Ua q)qy 



1 ^1 /-1-a 



dw / da„ 



Jo Jo 



da?; / da^w 



_pqX(a q ,aLl - a q - aL 



m c ~ (P + Ua ql) 2 



du [ 

II JO 



+ I du I da q J da' q 

Jo 



l-a' 



9 _ pqX(a' q , Og, l-a' q - a q ) 
? m 2 — (p + n(l — a g )g) 2 



(67) 



1/2 /-1/2 



JO 



If = e c (gg) / da q / da^ - - - ( e * - e ' ] S 7 (a g , a ? , 1 - a, - a q ) 



" 1/2 Z^ 2 2-a -3a- 
- < r (qq) I da q / da " ' 



jo (l-a g -ag) 2 v ' 



e * — e * i„, a„, 0:77, 1 — a„ — a 7 , 
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+ 2e c (qq) e « - e 



x 



1/2 /"^ , /" 1/2 1/2 -a- 

da q / da' / dag--- q —-rTHa q , a q , 1 — a q — a q ) 



+ / da, / ' da; / 9 da ? 1/2 a * T 



'■y V^9' ^ ^<j[ 



/•1/2 ]_ rag 

/ da 9 / da^— / da^Z y 4 (a g , ai, 1 - a q - oQ 
Jo Jo Za q Jo 



1/2 r a i r l ~ a 'q 

' I da i^f\ ^2 / da 'q / dag7^ 4 (a;, ag, 1 - a^ - ag) 

o A 1 — a q) Jo Jo 



(68) 



nl/2 nl/2 

?g = e q (qq) / da g / da^( e 
Jo Jo 



1 2 3q^^ Ota , \ 

-S + q - ^(T 3 - T 4 ) 



1 — a q — Ctq (1 — OL q — Ot^)' 



1 g Gt q 



Pi - r 2 ) 



a^ a^ 
(1 - a ? - a^) 2 



(a^, a^, 1 a^ 



+ 2e g (gg) e 



nil _£o\ — m 2 D , 
t — e * 1 



i-l/2 ra, i-l/2 1/2 -a- 

x W da 9 / da^ / da ¥ - q —^ [%_ - T 2 - T 3 + T 4 ] (a q , a q , 1 - a q - a q ) 

'O Jo Jo I 1 — a q ~ a q) 



1 nl/2 nl-a 



1/2 - aq 

1/2 JO ' Jl/2 {l-aq-aqf 



+ / da g / da^ j da q 



[Ti - T 2 - T 3 + %} (a q , a q , 1 — a q — a q ) 



nl/2 rl-aq | ra-g 



i -T 2 -T z + %} (a q , a' 1 - a q - aL) 
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1/2 

da c 



2(1 



a q ) 2 



da' 



l-a' 



dag [Ti - T 2 - T 3 + T 4 ] (a' a-g, 1 



(69) 
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